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We study elastic manifolds in a TV- dimensional random potential using functional RG. We extend 
to N > 1 our previous construction of a field theory renormalizable to two loops. For isotropic 
disorder with O(N) symmetry we obtain the fixed point and roughness exponent to next order in 
e = 4 — d, where d is the internal dimension of the manifold. Extrapolation to the directed polymer 
limit d — 1 allows some handle on the strong coupling phase of the equivalent iV- dimensional KPZ 
growth equation, and eventually suggests an upper critical dimension d u « 2.5. 



Disordered elastic systems are under extensive study 
both theoretically and experimentally. They are of in- 
terest for a number of physical systems, such as CDW 
j] , flux lattices 0, Q , wetting on disordered substrates 
4J, and magnetic interfaces [5|, where the interplay be- 
tween the internal order and the quenched disorder of 
the substrate produces pinned phases with non-trivial 
roughness and glassy features 6] . Typically they are de- 
scribed by elastic objects, with internal d-dimensional co- 
ordinate x, parameterized by a A^-component height, or 
displacement field u(x). Analytical methods are scarce, 
and developing a field-theoretical description poses a con- 
siderable challenge. One reason is that naive perturba- 
tive methods fail, technically due to the breakdown of 
the dimensional reduction phenomenon |7|, and physi- 
cally because describing the multiple energy minima in 
a glass seems to contain some non-perturbative features. 
One subset of these problems, the directed polymer (i.e. 
d = 1) in a random potential, maps onto the KPZ growth 
problem, well known to exhibit a strong coupling phase, 
which is out of reach of standard perturbative methods 
8]. It is thus important to obtain a field-theoretical de- 
scription of this phase, since the value and even the ex- 
istence of its upper critical dimension is still a matter of 
considerable debate H . 

One method which holds promise to tackle this class 
of problems is the functional renormalization group 
(FRG). Although it was introduced long ago, within a 
1-loop Wilson scheme [3 El 13, it is, not so sur- 
prisingly, hampered with difficulties, and only recently 
attempts have been made to push the method further 
El 111 IH El 13 HI IH Eg. The main problem is 
that the effective action at zero temperature becomes 
non-analytic at a finite scale, the Larkin scale, where 
metastability appears. Although fixed points are acces- 
sible in a d — 4 — e expansion, non-analyticity results in 
apparent ambiguities in the renormalized perturbation 
theory at T = [l9| • These problems are absent at 
T > [Hill (at least at leading order and for N = 1) 
but since temperature is dangerously irrelevant, the fi- 
nite temperature description is rather complicated |2l| . 
Until now, it has lead to a complete first-principle solu- 
tion of ambiguities (and calculation of the /3-function to 



four loop) only for the toy-model limit d = 0, N = 1 |22| . 
A case where ambiguities have been resolved from first 
principles at T = to 2-loop order, is the N = 1 depin- 
ning transition [lr! Il8j . Fi nally , the FRG has also been 
solved in the large- N limit [2fJ. Its solution reproduces, 
apparently with no ambiguity, the main results from the 
replica-symmetry-breaking saddle point of Ref. m > an d 
also underlies the importance of specifying the system 
preparation |2fj|. 

In the more difficult case of the statics within the 
d — 4 — e expansion, detailed analysis to two and three 
loops [11 13 |26j for the case of N = 1 have suggested 
several methods to construct a renormalizable field the- 
ory. These methods give a unique finite /3-function, with 
non-trivial anomalous terms. This /3-function satisfies 
the potentiality constraint, with anomalous terms dis- 
tinct from those at depinning, and a fixed point with the 
same linear cusp non-analyticity as to one loop, hence 
confirming the consistency of the picture. 

The aim of this paper is to extend these methods to 
the A^-component model. We show how an extended /3- 
function can be obtained and point out the specific fea- 
tures of the case N > 1. For the case of 0(A^)-symmetric 
disorder we compute the fixed point and roughness expo- 
nent £ to next order in e = 4 — d, where d is the internal 
dimension of the manifold. We then study the extrapo- 
lations to the directed polymer limit d = 1, and discuss 
the various scenarios for the strong coupling phase of the 
equivalent A^-dimensional KPZ growth equation. In one 
of them, a value for the upper critical dimension is esti- 
mated. 

We consider the model for an elastic Af-componcnt 
manifold 



H = J d d x^(Vu) 2 



V(x,u) 



(1) 



in a random 



potential with 
= S d (x - x')R{u - 



second cumulant 
u'), where u = u % 



V{x,u)V(x',u' 
is a A^-component vector. We derive general equations, 
and later focus on the O(N) isotropic case, noting 
R(u) = h[r) with r = \u\. Introducing replicas we obtain 



2 



the replicated action: 

^r = /d^E^'-^E^K-^) ( 2 ) 

a ab 

We now carry perturbation theory in the disorder and 
compute the one-loop and two-loop corrections to the ef- 
fective action T[u]. We use the usual power counting of 
the T = theory, identical to the case N = 1. Infrared 
divergences for d — 4 — e only occur in the 2-replica term, 
which at zero momentum defines the renormalized disor- 
der; there is no correction to the single replica term. The 
graphical rules are depicted in Fig. ^ We use functional 
diagrams, and mass rcgularization. The method and no- 
tations are identical to to which we refer for details. 
Here we only stress the differences with the case N = 1. 

The 1-loop correction to disorder (graphs a and (3 in 
Fig. 2) reads: 



S^iu) = Q[fiyi2(v)] a ~ d ij R(Q)d ij R(u) 



(3) 



Summation over repeated indices is implicit everywhere, 
and I = j k G\ = m- e I with G k = (k 2 + m 2 )" 1 . We de- 
fine the dimensionless function S 1 ^) := m e 5 1 R (recog- 
nizable by the parenthesis around the argument R). For 
later use we also denote the bilinear form 8 1 (R,R) := 
5 1 (R). This yields the standard 1-loop FRG equations, 
recalled below, and dijR develops a cusp non-analyticity 
at u — beyond the Larkin length scale L c . For the 
0(N) model one has d l] R = + inuj{h" - = 

h"(0)S l3 + \h'"{Q)r{5ij + u&j) + 0(r 2 ) and thus h"'(0) 
becomes non-zero at L c (u = u/\u\). 

The 2-loop corrections to disorder can be decomposed 
into a "normal" part, which is the complete result when 
R(u) is analytic |l5j . and an "anomalous" part which 
arises from non-analyticity. The normal part reads: 

S 2 n R(u) = (d l3 R{u) - d t3 R(0))d m R(u)d 3kl R(u)I A (4) 



-d l]U R{u){d, k R{u)-d %k R{Q)){d 3l R{u)-d iR(V)) 



I 1 



The first line stems from diagrams b and a of Fig. 1 re- 
spectively and the second from g,h,i,j. One has I a — 
SkiM G k 1 G k2 G 2 ki+k2 = mr 2t I A and we denote in analogy 
to S 1 (R) the dimensionless function <5( 2 )(i?) := m 2t S 2 R. 
The FRG /3-function is then: 

- md m R\ Ro = e[R + S\R) + 26 2 (R) - S^R)] , (5) 

where the repeated 1-loop counter-term S 1,1 (R) := 
2S 1 (R, 5 1 (R, R)) arises when reexpressing the bare dis- 
order Rq in in terms of the dimensionless renormal- 
ized one, defined as m € R, as detailed in [l^. From J2J) it 
reads: 

S^iR) = [(dijR-dijRiO^dijS^R) 

-dijS^R^odijR}! 2 (6) 
d l3 5\R) = d ijkl R(d kl R - d kl R(0)) + d lkl Rd ]kl R (7) 



• a _ Rgb 




FIG. 1: Graphical rules, one loop and two loop diagrammatic 
corrections to the disorder 



The property of renormalizability amounts to cancella- 
tion of the 1/e poles between the two last terms in (JjJJ 
using j = Ea an d I A - \p = N 2 (± + 0(e )) 16]. The 
/3-function jSJ) is obtained from J3J), Q and (JJ as: 



-md m R(u) = eR(u) + -[%i?(w)] 2 - d i3 R{Q)d i3 R{u) 



-(d l3 R{u) - dijR(0)) -d ikl R(u)dj k iR(u) - a i3 



(8) 



The cancellation works perfectly for the normal parts. 
Anomalous parts, to which we turn now, produce the 
last term. 

We start with the anomalous part of the repeated 
counter-term: 



where we denote the limits of small argument v - 

IMj ■= dikiR(v)djkiR(v)\ v ^o 

Vij ■= dijkiR(v)(d k iR(v) - d k[ R(0))\ v ^ 



(9) 



0: 



(10) 
(11) 



which, in general, are direction dependent. For a O(N) 
model, the third derivative tensor: 

dij k R{v) = A(r)(5ijVk + 8ikVj + 6kjVi) + B(r)viVjVk (12) 

with v = v/\v\, A(r) = (rh"-h')/r 2 andB(r) = {r 2 h"' - 
3rh" + 3h')/r 2 , has a -C-dcpendcnt small v limit (|12() with 
A(0) = -B(0) = h"'(0)/2. This yields: 



fJtij 



N+ 1 



(13) 



and, similarly one finds z/y = ^^-h'"(0) 2 (5ij — ViVj). 

Let us first superficially examine the structure of the 
2-loop graphs, following the discussion in [J^- As for 
N = 1, one can discard c = d = from parity and 
similarly set m + n = and p + q = 0. One can then 
write: 

8^ (R) = -{fiah + h 3 I 2 )d l3 R{u) (14) 
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where the first term comes from graphs e (more properly, 
from the sum of all graphs a to /) and the second from 
graphs k + I (from the sum of graphs i to I). Global can- 
cellation of the 1/e pole in the (3- function works provided 
jlij + 2&ij = fly + Vij. This then produces = fiij/2 in 
the FRG equation above. 

We can now use the methods introduced in to an- 
alyze the total 2-loop contribution to the effective action, 
including possibly ambiguous graphs. One first computes 
T[u] in a region of u where no ambiguity is present, us- 
ing excluded replica sums, and constraints valid in the 
zero-temperature theory (th e so-called sloop elimination 
method, Section V.B in 19]). One finds that extraction 
of the 2-replica part yields /2,j = /i^ and 29ij = My, 
i.e. it works as for N = 1. This is equivalent to renor- 
malizability diagram by diagram, and thus it satisfies 
the global renormalizability condition. The background 
method also yields that result (0], Section V.C). The 
end result for the /3-function, oeij = fUj, although unam- 
biguous for A = 1, needs further specification for TV > 1, 
since the limit in (|13fl is direction dependent. 

Another important consideration for the resulting (3- 
function is the issue of the "super-cusp". For N = 1 it 
was found that the /3-function is such that the cusp non- 
analyticity of R"(u) at u = does not become worse at 
two loops. That by itself constraints the amplitude of the 
anomalous term, since any other choice yields a stronger 
singularity [2Sj. We now point out that if v and u, in JjJ}, 
(JTJ3J, (T2J, are colinear, i.e. fiij(v) — Hij(u) then there is 
no super-cusp. Indeed the result: 

Oij(u) = lim ^d ik iR(ru)dj k iR(ru) (15) 

obviously yields cancellation of the linear term in u in (|5|) 
(although it is not the only possibility H^J). Colinearity 
of v = u a — u c and u = u a — Ub is natural if one com- 
putes the effective action in a background configuration 
breaking the rotational symmetry, which appears to be 
required for the present theory to hold. 

We now specialize to the 0(A) model. Starting from 
(|5|) and further rescaling h(r) — > m~ 4 ^h(rrn^), using £ 
we obtain the following FRG flow-equation to two loops: 

-md m h(r) = (e - 4()h(r) + (rh'(r) 
h"(0)h"(r) 



+i(h"(r)-h"(0)) h'"(rf 
N-l (h'(r)-rh"(r)f (2h'{r)+r(h"(r)-3h"{0))) 



2 

-h'"{0) 2 



N - 



Mm a- 



1 h'(r) 



(16) 



where the last line arises from the anomalous term (I15JI . 
This FRG equation admits for any N a non-trivial attrac- 
tive fixed point such that h"{r) has a linear cusp at the 
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FIG. 2: Numerical results for the exponents C,\ and £2 for 
different values of JV (top). Numerical plots of £i(N) (bot- 
tom/left) and & (N) (bottom/right), in blue with the numer- 
ical values from the table as dots. The red curves (no dots) 
represent the asymptotic expansion. 



origin and decays to at infinity faster than a power law, 
thus corresponding to short range (SR) disorder. Find- 
ing the associated C — Ci e + C2^ 2 + 0(e 3 ) is an eigenvalue 
problem, which has to be solved order by order in e fol- 
lowing |la. llSl [jjj . Our results for d and £2 are given on 
Fig. |21 Although for SR disorder no analytical expres- 
sion can be found for £1 and £2, their large- TV behavior 
can be obtained from an asymptotic analysis of p6|l . Let 
us extend the analysis of Balents and Fisher (BF) [T^ |. 
Define h = 1/Nh, y = r 2 /2 and h(r) = Q(y). For y > 1 
the FRG equation can be linearized: 

(e - 2()Q' + 2QyQ" -{A + W)Q" - 2ByQ"' = (17) 



with A = (1 - ±)Q> + ^h"'(0) 2 and B = ±Q' Q + 

^h'"(0) 2 , h"{0) = Q'(0) = Q' . BF noted that there 
is an overlapping region 1 3> y S> N where the solution 
can also be found perturbatively by expansion in 1/N, 
yielding for Q a pure exponential. It is indeed an exact 
solution of (|17fl . with a unique value for £1, the BF result 
Ci « Ci with Ci = 1/(4 + AO (i-e. the result from the 
replica variational method [25j). The corrections (which 
arise from the neglected non-linear terms) are shown to 
be exponentially small; a more accurate estimate being 
Ci w Ci 1 with d 1 = Ci° + (A + 2) 2 /(A + 4) 2 2-( A, + 2 )/ 2 /(4e). 
To next order we find similarly the approximation to £2 



(A 2 - 1)(2 + A) 
2(4 + A) 3 (3 + A) ' 



(18) 
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where we have not attempted to estimate further correc- 
tions, presumably again exponentially small at large N. 
We note that C2 

arises from the anomalous terms only. 
These estimates are listed and plotted on Fig. together 
with the numerical solution of (|16fl . The quality of the 
large- N analysis is quite remarkable. 

We now discuss the extrapolation of our result to the 
directed polymer (DP) case d = 1, e = 3, plotted in 
Fig. |3 We see that the 2-loop corrections are rather big 
at large N, so extrapolation down to e = 3 is difficult. 
However both 1- and 2-loop results as well as the Pade- 
(1,1) reproduce well the two known points on the curve: 
C = 2/3 for N = 1 8] and C = for N = 00 20]. This 
branch in Fig. [3| corresponds to zero temperature and a 
continuum model. On the other hand we find that for all 
curves in figure [3] the roughness £ becomes smaller than 
the thermal Cth = \ at N — N uc ss 2.5. This naturally 
suggests the scenario that at non-zero temperature £ = 
1/2 for N > iV uc , i.e. N uc is the upper critical dimension 
[2( . The same argument gives an upper critical dimension 
N U g for the KPZ-equation of non-linear surface growth 
[H, [ll|. On the other hand, simulations on discretized 
models of both the directed polymer (at T = 0) and 
the KPZ equation 0, suggest that C > 1/2 in all 
dimensions, but should be taken with caution 30]. Since 
the FRG is a systematic expansion in e = 4 — d, such a 
scenario seems reconcilable with our above results only 
through non-perturbative corrections in e, possibly non- 
analytic at e = 2. 

To conclude we have obtained for the iV-component 
model a FRG description at 2-loop order. Various stud- 
ies, including at large N, are under way to obtain a bet- 
ter understanding of the structure of the theory. For the 
KPZ growth and the directed polymer we have improved 
the determination of the possible upper critical dimen- 
sion. Further numerics, in particular for the directed 
polymer at T = would be helpful. 
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FIG. 3: Results for the roughness £ at 1- and 2-loop order, as 
a function of the number of components N. We both show a 
direct extrapolation and the Pade (1,1): Cpa.' — — — 
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